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Abstract 



1 Introduction 

1.1 Summary 

The goal of this article is to expand the constructions presented in [1] for the 
case of hyperelliptic curves (case N = 2) to the general spectral curves of any 
order. 

As in [1] our main tool is the affine model of the Jacobi variety of the 
curve or, in the other words, the construction of the separated variables for 
the underlying integrable system. Although the existence of such construction 
has already been proven in [3] it only becomes really useful in the study of the 
integrable systems if it can be specified explicitly. It has been known for a long 
time for the hyperelliptic case (see [4] for example) and was originally done by 
Sklyanin in [5] for SX(3). His ideas based on the functional Bethe Ansatz were 
later generalized to SL(N) in [6] while a more elegant algebraic-geometrical 
approach (which is the closest to our one) was used in [7]. 

However neither of these models could be used to generalize the results of 
[1] because we need to choose a particular orbit carefully (see the section 3) 
and only then the construction of the separated variable gives us the additional 
insights into the structure of the cohomology groups. 

°Membre du CNRS 
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With this problem behind us, the other results obtained for hyperelliptic 
curves carry over to the general case with only technical problems. But in 
other respects the general case is quite different and almost surely much more 
complicated. Because of this, we are still unable to even conjecture the explicit 
form of cohomologies for it. We do show that the highest nontrivial cohomology 
group doesn't suffice any more to generate the elements of the ring of observables 
A for N > 2 in a marked difference with the hyperelliptic case. 

1.2 The Plan 

In the first section after this one we introduce the spectral curve £ which is the 
starting point of our study. In the context of integrable systems, this curve is 
the phase space of some integrable model and the the ring A of meromorphic 
functions on the Jacobian of this curve with singularities only on the theta 
divisor corresponds to the ring of observables. 

In the section 3 we construct the matrix space which will be later shown 
to provide an algebraic model of A and supply it with the Poisson structure. 
This task is quite difficult in the general case as it is less obvious to identify the 
'constant' ring elements (in the hyperelliptic case we only had to deal with the 
center of mass motion one) which must not be taken into account if we are to 
get the isomorphism with the observables ring. 

Next, in the section 4, we construct the separated variables which provide 
us with the mapping from the matrices lying in the space discussed above to 
the point on the Jacobian of E. 

Section 5 deals with a much more difficult problem of constructing the inverse 
mapping. Although it is not as explicit as in hyperelliptic case, we still were 
able to find a constructive way to express it. The existence of this mapping 
allows us to conclude that our space of matrices is indeed isomorphic to the 
affine Jacobian. 

Finally, we combine the preceding results together in the section 6 to cal- 
culate the Euler characteristic of the Jacobian and discuss the cohomologies 
structure which would be compatible with it. As indicated above, we are unfor- 
tunately forced to conclude that the simplest situation encountered in studying 
the hyperelliptic case doesn't hold in general. 

2 The Spectral Curve 
2.1 General properties 

We consider an algebraic curve £ of genus g and its theta divisor X which is a 
(g — 1 )-dimensional subvariety of the Jacobi variety J associated to £. We call 
the affine Jacobian of E the non-compact variety J — X. 

Our construction of the algebraic model of the affine Jacobian is only appli- 
cable to the spectral curves, that is curves defined by the equation 

r(w, z)=w N + h {z)w N - 1 + . . . + t N (z) = (1) 

where tk(z) are polynomials in z of degree at most kn — 1 and the degree of 
t]y(z) is exactly Nn — 1. 
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The genus of the curve (1) is easily computed to be 
9=\{N-l)(Nn-2) 

and we can write down the abelian differentials of the first kind on it explicitly 

as 

, > w l z k , . , f < l < N - 1 

akl{w > z) = d^7) dz Wlth \0<*<(JV-/-l)n-l 

(we shall also denote them below as 07 with a composite index / = l,...,g 
instead of two indices k and I for convenience) . 

Also note that the second kind differentials on £ have the same form except 
that I and k are unlimited for them. 

Let us also choose a canonical homology basis of £ (ot\, . . . , a g , (3i, . . . , (3 g ) 
and define the corresponding normalized differentials ujj 

9 
J=i 

where A is the matrix of a-periods: 



i/j = f <?i 

I a j 

These differentials are normalized because we have 

f § ai uj - 5 U 

where B is the period matrix (it is symmetric, with positive definite imaginary 
part). 

The Jacobi variety J of S is then simply defined as a g-dimensional complex 
torus 

U> + BZ9 

The Riemann theta function associated with the Jacobian J is defined by 
°(0 = exp2Tri t n(^Bn + 
It satisfies the quasi-periodicity condition 



0(£ + m + Bn) = exp 



-2iri t n(^Bn + ^) 



0(0 



for any m, n G Z 9 . 

Thanks to the above we see that the 0-divisor 



© = {mo = 0} 

is a (possibly singular) (g — l)-dimensional subvariety of J. 
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2.2 Holomorphic functions on the affine Jacobian 

From the point of view of the integrable systems theory, the observables should 
be the holomorphic functions on the affine Jacobian or, in other words, the 
meromorphic functions on J with singularities on the theta divisor only. 

Such functions may be described using the products and quotients of the 
theta functions themselves but for our purposes an alternative description is 
more convenient. 

First, let us introduce the Abel map. Note that the curve (1) has a single 
point at infinity which we denote by oo as in local coordinates z' — 1/z and 
w' = w/z N the curve equation takes form w' N + z'(. . .) = 0. Accordingly we 
may define the Abel map using oo as the base point: 

A:Y, x a _> J 

(pi, ■•■,?,) - ELiC^ + A M 

where A is the Riemann characteristic (a vector of constants depending only on 
E and the choice of the base point). The map A is well-defined as deforming 
the integration path by a loop ai only changes the fc-th component of the RHS 
by Ski while making an extra (3i loop adds Bki to it and thus deforming the path 
in any way doesn't change the point on the Jacobian. 

By a corollary of Ricmann's theorem any (e6 may be written as 

r q 3 

i=l J °° 

If we write £ as an image of a point of S Xff by the Abel map: 

£ = A(p U ...,Pg) 

we see that we have two equivalent divisors: 

(pi,...,p g ) ^ (qi,...,q g -i,oo) 

or, in other words, there exists a function with the simple poles in pi and simple 
zeroes in qj and oo (we can always shift it by the value in oo to ensure that it 
has a zero there). 

If all pi are distinct, to find a function with the singularities on 6 we have 
to find a function with simple zeroes in (pi). This is possible if and only if 

A(pi, ...,p g )= dot (wi(pj)) = 

The situation is more complicated in the case when two (or more) points 
coincide but we can cancel the zeroes due to this by considering functions of the 
following form: 

det (a(pj)) 

A(pi,...,p fl ) {> 

where a are any differentials of first or second kind. 

Thus in any case the functions of the form (2) have poles on the theta 
divisors only. Moreover, by reversing the argument above, we also see that we 
can get all holomorphic functions on the affine Jacobian in the form of products 
of functions (2). 
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3 The Affine Model 

3.1 The Matrix Poisson Brackets 

Consider a space of N x N matrices of polynomials of z: 

M(z) = {mijiz))^^ 
We start with the usual expression for the Poisson brackets for such matrices: 
{m(zi) ® m(z 2 )} = [r 12 (z 1 ,z 2 ),m(z 1 ) + m(z 2 )] (3) 
where the r-matrix is defined by 

1 fl 

ri2{z u z 2 ) = - {zi + z 2 )t oa + Zl r + + z 2 t + 

z\ — z 2 \Z 

and 

i 00 = ^EuQEii 

t + - = ^ E V ® E Ji 

i<j 

t~+ = E V ® E H 

i>j 

where the matrices E^ are N x N with the only non zero element being 1 at 
the intersection of i-th row and j-th column. 

The brackets (3) are compatible with the following form of the polynomials 

m(z): 

j m,ij(z) = wbf^ z n ~ x + . . . + m[™~^ if i < j 
1 mij(z) — mnf^ z n + . . . + m£™ ^ z otherwise 

With such form of the matrix elements we can easily see that 

detm(s) = z Nn - 1 ^ JJmf^j + 0(z Nn - 2 ) (5) 

because only the minor containing the last element of the first row and diagonal 
below the main one contributes the (Nn — l)-th power of z. 

The characteristic polynomial of m(z) has the following form: 

det (m(z) + w) = w N + t 1 (z)w N - 1 + ... + t N (z) (6) 

with coefficients tk(z) being polynomials of degree kn — 1 in z: 

t k{z ) = z ^-HW + ... + ti kn -v 

which means that the equation 

det (m(z) + w) = (7) 



(4) 
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defines a spectral curve in the sense of the previous section. We shall fix the 
characteristic polynomial of m(z) so that this equation describes the same curve 
£ as the equation (1). 

All coefficients are in involution which follows from the identity 

{logdet(w + m(zi)),m(z 2 )} = 

= tY®l([{w + m(z l )®I)y 1 [r 12 {z 1 ,z 2 )J®m{z 2 )]} (8) 

(which is easily proved by replacing log det with tr log and developing the brack- 
ets of the logarithm). 

Moreover, many of these coefficients lie in center of the Poisson algebra: 
indeed, the coefficients . . . tjj." ^ appear in front of degrees of z which are 
too big to appear in the RHS of (8). In particular, the dominant coefficient of 
the determinant (5) lies in the center. 

3.2 The 'Physical' Matrix l(z) 

We cannot hope that m(z) provide us with an algebraic model of the affinc 
Jacobian for several reasons. 

First and the most obvious one is that we want to embed a curve of genus 
g = |(7V— l)(A?n— 2) into the matrix space whose dimension is, as seen from (4), 
equal to nN 2 . Choosing the subspace of matrices with the fixed characteristic 
polynomial yields n + 2n + . . . + Nn = ^nN(N + 1) constraints somewhat 
reducing the dimension of this subspace but we still have 

nN 2 - ^nN(N + 1) = ^nN(N - 1) = g + (N - 1) 

and so we must get rid of N — 1 extra 'unphysical' degrees of freedom if we are 
to establish an isomorphism. 

There is another motivation for doing what follows which is less obvious but 
even more important: in fact, it is impossible to reconstruct the full matrix 
771(2) in terms of the single- valued functions on the jacobian as it contains some 
^-function like pseudo functions which cannot be expressed with the help of 
holomorphic functions on (a subvariety of) J . It also explains why eliminating 
these extra degrees of freedom is far from being trivial: their commutation 
relations with the other elements are quite involved (in particular, they are not 
at all in the center of the Poisson algebra) and so we cannot just remove them 
from consideration. Instead we have to go through the construction below in 
order to obtain something that we can expect (and actually will) to express in 
terms of functions on the jacobian. 

To solve both of these problems, let us replace m(z) with a conjugated matrix 

l{z) — sm(z)s^ 1 

with s built from m(z) in the way described below: 
We have, because of (4), 

m{z) = z n ii- +z n - 1 ii + 0{z n - 2 ) when z 00 (9) 
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with fi~ being a strictly lower triangular matrix and fi - an upper triangular 
one (with diagonal terms). Then if we denote the first row of fi as fi\, we define 
the N x N matrix s as 



V 



mi 



(note that it is a lower triangular matrix - simply because pT is) . 

Replacing m(z) with l(z) plays the role of 'gauge fixing': indeed, l(z) still 
defines the same curve (7) so the physical picture doesn't change but all the 
extra degrees of freedom have disappeared from l(z). To see it, let us look at 
the degrees of the polynomials l(z)ij. 

First, we have: 

l{z) = s f i- S - 1 z n + 0{z n - 1 ) 
= Uz n + 0(z n - 1 ) 



where 



U 



( o 
1 




o \ 





(10) 



V o 10/ 

which shows that the degrees of all elements except the ones on the diagonal 
below the main one are strictly less than n. 

Next, let us determine the degrees of the elements in the first row. As 



we have 



Iin 



{sfis %■ 



0(z n - 2 ) 

V JV- 



sn(M s_1 )ij 
sii(ss _1 ) W j 
5NjSn 



U<N- 1) 



= Sll z- - 1 + 0(z n - 2 ) 
So the degrees of the elements of l(z) have the following form: 



/ n- 2 n- 2 n- 2 
n n — 1 n — 1 
n — 1 n n — 1 



deg z /(z) = 



\ n-l 



n-2 n-l \ 
n — 1 n — 1 
n — 1 n — 1 



(11) 



n - 1 / 



Note that sn, which is the dominant coefficient of Zijv ; is also exactly the 
dominant coefficient of the (5) and so it lies in the center of the Poisson algebra 
or, in other words, is a constant. 
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Finally, the dominant coefficients of the other elements of l(z) with maximal 
degrees, that is the polynomials hi(z), ■ . . , In n-i(z) are constants as well (as 
seen from (10). 

Calculating the number of coefficients in the matrix l{z) using (11) (not 
counting the constant ones determined above) we see that now we have exactly 
the 'correct' number of parameters. 

3.3 Calculation of the Poisson Brackets of l(z) 

Our next goal is to calculate the Poisson brackets of l(z). For this we start by 
computing the brackets of fj, and m(z). 

First, fix z 2 and let z\ — > oo. Then we have: 



{/if m(z 2 )} 




(12) 



with 



Next step is to develop 




(13) 



where 



«12 = (I <S> s) {s f m(z 2 )} (s" 1 ® s _1 ) 



To compute k\ 2 we first write 



{ Sj fm(z 2 )} = { Ml ®m{z 2 )} (( A1 -) Af +( m ®/){( M -) Ar j ®m(z 2 )} 



where Sj = \i\ [p, ) 



N-j 



is the j-th row of s. Using (13) we have 



{ Sj ® m(z 2 )} = 
= ® £ n ) ((m - )^ ® /) + *2 E i k 0"~) 



JV-j+l 



(Mi®/)(^ 00 +*- + )((M-) JV " j ®/) + 

+ (^i®/)(^ 00 +t- + )((M-) JV " J '®/)- 
( Ml ®/)(( A1 -) Ar - J ®/)(>+t-+) 



Hence after simplification 



{s ® m(z 2 )} = (s ® 7)(--i 00 - r + + -I ® E n ) + z 2 £ £7 jfc (jO 



JV-j+l 



(giEfei 
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And so 

K12 = -(s®s)(it 00 +t~+ -^ItoEu^s- 1 ®s~ 1 ) + 

z 2 {I®s) [^E ok {^-) N -^ l ®E k ^j (s" 1 ®*- 1 ) 

We can leave out the term with / ® En as it commutes with l(z\) ® I. The 
second term can be simplified using the multiplicative properties of E rs : 



(l®s)[Y,E lk (P~) N ~ 3+1 ®E kl \ = J2E jk (ti-) N - j+1 ®s lk E n 

jk J jkl 



ji 

Y^E ol U N+1 -i ®E^j (s® J) 

with the matrix U defined in (10). 

Finally, inject this expression into (13) to obtain the Poisson brackets of l(z): 

{l( Zl ) ® l(z 2 )} = [fi 2 («i, z 2 ), l(zi) ® 7] + [f 21 (z 2 , zi), I ® Z(z 2 )] (14) 

where __ 

fi2(^i,%) = ^^t + z 2 J2E jk U N+1 - j ®E kl 
zi — z 2 ^— ' 

The equation (14) specifies the Poisson structure on the matrix elements 
which we are going to use henceforth. 

Note that we could have started directly from this expression of the Poisson 
brackets, however obtaining it from the considerations used above gives some 
additional insight into why do we have to use the r-matrix f\ 2 here. 



4 Separation of Variables 

In this section we are going to construct the separated variables on the Jacobian 
from the matrix l(z). This is the first step of our plan to establish an isomor- 
phism between the (affine) Jacobian and the space of these matrices, although 
it also provides a generalization of the construction of the separated variables 
originally proposed in [5]. 



4.1 The Baker— Akhiezer function 

Consider the function tp which is the eigenvector of l{z) for the eigenvalue — w: 

{l(z) + w)i>(w,z) =0 (15) 

Any of the components of the vector tp( w i z ) is a Baker-Akhiezer function. 
The poles of ipi(w,z) are fixed (they are the points of the dynamic divisor 
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describing the time evolution of the system at the moment t = 0) and the 
zeroes of them give the separated coordinates as we'll see below. 

Let us choose the first component and let (wi,Zi) be the g zeroes of ipi(w 7 z). 
Then write l(z) as: 



l(z) 



' hi(z) 


b(z) 


* 






d(z) 


V * 





with the notations inspired by N — 2 case where the matrix was written as 

a(z) b(z) 



l(z) = 



c(z) d{z) 



Let us also split tp likewise: 



ip(w,z) 



( ipi(w,z) ^ 
ip 2 (w,z) 



4>{w,z) 



V 1 / 

Then for any Zi decomposing the equation (15) defining tp gives 

( b(zi)<f)(wi,Zi) = 

\ d(zi)(f>(wi,Zi) = -Wi4>(zi) 

Or, combining them: 

b(z i )d k {z i )<j){w i ,Zi)=Q Vfc>0 

This leads us to the idea to consider the matrix 

/ b(z) \ 
b(z)d(z) 



(16) 



Z(z) 



V b(z)d N ~ 2 {z) J 

Zi) — for any i and hence, as (f> ^ 



(16) implies, of course, that Z(zi)4>(w. 
(its last component is 1) we have 

detZ( Zl ) = \/i = l,...,g 

Remembering the specific form of degrees of the matrix elements from (11) 
we can see that the maximal degree of z in det Z(z) comes from the product 
containing Zuv and the elements from the second main diagonal in the matrix 
l(z). This degree is equal to J2k=i {kn — 1) = g and so the determinant det Z(z) 
is a polynomial in z of degree g which has g roots zi which proves that it is equal 
(up to a constant factor) to Yl^z — Zj). 

This allows us to now define Zi as the roots of the polynomial det Z(z) 
which is constructed from the matrix elements alone (and without any further 
reference to the Baker- Akhiezer function). Again, we could have started by 
defining them in this way as the subsequent calculations only use the properties 
of the matrix Z but making the link with the Baker- Akhiezer functions makes 
it more clear why do we need to define Zi in this way. 



10 



4.2 Determination of m, 



Next we are going to show that Wi are unambiguously determined by the matrix 
elements as well. 

For this, let us introduce the N x (N — 1) matrix 

b(z) 
d{z) + w 

Considering its action on the vector 4>(iVi, Zj) we see that the N determinants 
obtained by removing the fc-th row from this matrix for k — 1 , . . . , N are null. 
Thus, for any (TV — 1) x A*" matrix A we have 



det A 



b{ Zi ) 
d(zi) + w t 



= 



(17) 



Let us apply this with 



/ 1 0...0 

b(zi) 

b{ Zi )d{zi) 

b(zi)d N - 4 ( Zi ) 



where £ is an arbitrary row vector. 

Then (17) yields, after eliminating the other terms by line subtraction, 



det 



/ Kzi) \ 

b(zi)d(zi) 

b( Zl )d N - 3 (z t ) 
V t(d{zi)+Wi) J 



(18) 



At this point we need to make an additional genericity hypothesis: let us 
suppose that Zi is a simple root of det Z(z). This ensures that the matrix 

( b{zi) \ 
b(zi)d(zi) 

\ b(Zi)d N - 3 (Zi) ) 

is of rank N — 2. Indeed, if this were not the case we would have had 

JV-2 



J2Hz l )d k - 1 (z l ) = 



k=i 



and hence 



JV-2 



J2b(z)d k -\z) = 0(z-z i ) 



k=i 
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and also 



N-2 



fe=i 

but the last two equations imply 

det Z(z) = 0((z- Zl f) 

which contradicts our hypothesis. 

Thus this hypothesis ensures that there exists a £ for which wi can be defined 

by 

/ Kzi) \ 
b(zi)d(zi) 



det 



W, = -- 



b(Zi)d N - 3 (Zi) 

td(zi) 



( b( Zi ) \ 
b(zi)d(zi) 



(19) 



det 



b(Zi)d N - 3 (Zi) 

V e / 

as, clearly, then (18) is satisfied. 

So we have achieved our goal of constructing g points (wi, z{) from the given 
matrix elements. 

Note that although Wi can still be uniquely constructed even in presence 
of double roots of det Z(z), we would have other, more serious, problems with 
reconstructing the matrix elements in the next section without the genericity 
hypothesis and so we will assume it from now on. 



4.3 The Poisson brackets 

So far we have seen that, starting from the matrix elements, we may define a 
point (v)i,Zi) on the jacobian. But it still remains to show that w and z defined in 
this way satisfy the canonical commutation relations and so are really separated 
variables in our problem. 

These commutation relations have been first calculated in [7] but, as men- 
tioned in the introduction, we propose a more algebraic way of computing them 
which allows a relatively straightforward generalization to the quantum case as 
seen in [8]. 

So in this section we are going to carry out the calculations proving that 
we indeed have the expected commutation relations for the variables we have 
just constructed. To simplify them, first let us note that we can define Z(z) in 
terms of the matrix m(z) instead of l(z) as conjugation by s doesn't change the 
brackets of the determinants but the calculations which follow become much 
simpler because of a simpler form of the brackets of m(z) itself. 

Thus in the rest of this section we are going to redefine b(z) and d(z) as the 
corresponding parts of the matrix m(z) and not l(z), in other words we have 

/ m n (z) b(z) \ 

m(z) = 

V : d(z) J 
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where b(z) a covector of size N — 1 and d{z) is a square matrix of order N — 1. 

To abbreviate the notations, we will omit the argument z\ and use primes 
for the functions of z 2 . We will also use the short hand notations: 

x\ = x(z\)®I — x®I 
x 2 = I®x{z 2 ) = I®x' 

and 

{x{ Zl ) ® x(z 2 )} ijM = {x U X 2 }ij,kl = {Xik,x'ji} 

for any x. 

Finally, we introduce a special notation for the rows of the matrix Z(z): 

a ( - k) {z) = b{z)d k - 1 {z) 

4.3.1 The brackets {z i} zj} 

Our first goal is to show that 

{dct Z, dctZ'} = 

For this first let us use the identity 

{dct Z, dct Z'} = dct Z dct Z' (tr ® tr) (Z^ 1 Z 2 1 {Z 1 , Z 2 }) (20) 

which shows that we need to calculate {Z±, Z 2 } and for this, in turn, we need 
to know {61,62}, {bi,d 2 } and {di,d 2 }. To compute them we start with the 
definition (3) 

{mi,m 2 } = [r(z 1 ,z 2 ),Am(zi,z 2 )] (21) 

where 

Am(zi, z 2 ) =m\ +m 2 = m®I + I®m' 

and split this relation according to the decomposition of the matrix m(z): let 
us introduce the matrix p which is the same as r(z\, z 2 ) but of order N — 1 (to 
be precise, Pij,kl = r i+i j+i,k+i i+i)- Then we immediately get from (21): 

{di,d 2 } = \p(z 1 ,z 2 ), d\ + d 2 ] 

As Am(zi, z 2 )ijM = unless i — k or j — I and r^^i ^ only if i = I and 
j = k we see that (21) also gives 

{61,62} = 

And using the explicit form of r we obtain 

{buck} = — — — $6(2:2) - hp 
zi — z 2 



The object $ appearing here is a notational convenience and is formally defined 
as 

N 

$ = ^ e i ® fi 
i=l 

with ei and fi being a row and column matrices respectively with 1 at i-th 
position and elsewhere. $ transposes the matrix dimension in the tensor 
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product, i.e. when it is applied to b 2 which is a matrix in the first component 
and a row in the second one it yields a tensor product of a row in the first 
component and a matrix in the second one defined by 

{®b 2 )i,jk = Mfc 

We also define '$ which results in a tensor product of a matrix and a row and 
which acts on b\ as 

(*$&i)i,jk = Sikbj 
Our next goal is to calculate {b\d k , b 2 d l 2 } for which we need 

(as the other terms simply cancel pairwise) and also 

{b u d i 2 } = -^drv^+E^^^M^ 

p— 1 p— 1 

/c k 
{dib 2 } = -b 2 Y / dr l pd h r q +Y / dr i ^^ t <s>b 1 d k 1 - q 

g=l g=l Zl ~ Z2 

Let us define some more abbreviations: 



d[f = E k g zUr i P d k 1 - q - 1 



and 



Putting all the above together and using the new notations we obtain: 

{fcidf- 1 ,^ -1 } = 
= b 2 {&!, 4" 1 } d*- 1 + 6x {d*"\ 6 2 } 4" 1 + 6x6 2 {d k ~\ 4" 1 } 

+ y 4 p) _ ^—^'-^df- 1 + y ^— «#a< fc - 9) 4- 1 (22) 

Now we have all ingredients to calculate the brackets of Z(z) itself: for this 
we just rewrite (22) in a matrix form: 

{Zi, Z 2 } = —Z\Z 2 \ — Z 2 Z\ 2 + 

+ —^—y{z 2 Z' rs {D s ®U r Z') + z 1 Z rs {U r Z®D l s )) (23) 
Zi — z 2 — ' 
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where Z 2 \ is a matrix whose rows (relative to the second tensor product com- 
ponent) are a^: 



JV-l 



£21(21,22) = X! i 1 ® fk)a. 



(k) 
21 



fc=i 



and £12(21, 22) is defined in the same way using a( 2 '■ 

JV-l 



£12(21,23) = ^ (/; <g> J) a 



Finally, D r is just the matrix tailored to allow us to pick the r-th element 
of d in the matrix form: 



D r {z) = 



( & r \ 

e r d{z) 
\ e r d N - 2 (z) ) 



and the usage of U (defined by (10)) lets us to sum over all values of r instead 
of just the ones in range 1 < r < i (for matrix element ij, kl). 

We can now inject (23) into (20). Let us consider the four terms resulting 
of (23) independently, starting with the first one (containing Z 2 i) which results 
in: 

det Z det Z' (tr <g> tr) {z 2 x Z 2 \{z 2 , 21)) 

To calculate trig)/ note that only p contributes and hence, using its explicit 
form: 

(ir®I)p(z 1 ,z 2 )= l - Z -^^I 

I Z\ — z 2 

so, leaving the factors aside, we have only 

{I®tv)Z^ l Z 2l 

to calculate. And to do this, we simply notice that, after eliminating p from 
d 2 k i we are left with (k — l)b 2 d( k ~^ or, in other words, the (k — l)-th line of Z 2 
times k — 1 and so this trace reduces to 

tr (KZ 2 Z^) 

with 



K 



( 
10 
2 

Vo 



... \ 




TV -2 J 



(24) 



which of course means, as tr K = 0, that the contribution from the first term is 
0. The contribution of the second term vanishes in exactly the same way. 
Now consider the third term which contains a sum of 



Z' rs {D s ® U r Z') 
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The only contribution to the trace over the second space, after Z' cancel, is 
tr U r which yields a as a factor and so this term vanishes as well - just as the 
fourth one does for the symmetrical reasons. 
So we have finally established that 

{det Z, dot Z'} = 

which proves that Zi defined as the roots of det Z(z) also satisfy 

{zuZj} = Vi, j 



4.3.2 The brackets {zi,Wj} 

The next step is to calculate the (only non trivial) brackets of z and w. Let us 
define B{z) and A(z) by 



B( Zi ) = 

Wi = —A(zi) 



i.e. B(z) = det Z(z) 



As we have seen above, B(z) is a polynomial in the matrix elements and A{z) 
is a rational fraction in them (defined by (19)). We use these notations just to 
emphasize that these B and A are the generalizations of the functions from [5] . 
Then we have 



and so 



0={B(zi),Wj} = {B(z), Wj }\ z=Zi +B'(z i ){z i ,w j } 
{B(z),A( Zj )}\ z=Zi 



{Zi,Wj} 



B'(zi) 



(25) 



which means that our task is to calculate {B(z), A(z')}. 

For this we will reuse the results of the previous subsection: note that A(z) 
can be written as 

., , detX(z) 
A(z) - 



with 



X 



( a« 

a (N-2) 

V Crf J 



det Y{z) 



and Y = 



a 



(JV-2) 



V s 7 



and so to calculate the brackets in question we have to calculate the brackets of 
Z with X and Y which is done similarly to the calculation above with the only 
differences being in the last row. 

Accordingly we are going to proceed in the same order as above so as to 
obtain a form similar to (23). Thus we first define the matrices corresponding 
to Z: 



-^21 = 



; (JV-2) 



and Y21 = 



V bp ) 



-AN-2) 

hi 
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Note that the rows here are relative to the second tensor product component, 
&2i already being matrices in the first space. And X\ 2 and Y\ 2 are defined in 
the same way using d\ 2 and £1 for the last row (relative to the first space) of 

X\2. 

Let us start with the easier computation, that of {Z\,Y 2 }: in this case we 
already know the result as it will be the same as (23) for all rows (relative to 
the second space) except for the last one for which it is simply null. So we can 
write 

{Z U Y 2 } = -Z 1 Y 21 -Y 2 Z 12 + 

+ ( Z2Y rs ( D > ® ^ Y ') + Z ^ Z rs (^'Z ® ID')) ) + 



+ (/®/jv-i)6^i 



(26) 



with I and U being, respectively, the identity matrix and U without the last 
row: 

/ ... \ 



I = 



( 1 



V 



\ 



u = 



1 
1 



V o 






J 



Note that U r Z = U r Y = U r X for any r anyhow and we only use U r Y above 
for cosmetic reasons. 

Using again the identity (20) we get only the contribution from the last term 
of (26): 

{det Z u det Y 2 } = det Z x det r 2 tr 2 (1 <g> /jv-i) V^ 1 ( 27 ) 

with 

as all the other terms vanish inside tr <g> tr for the same reasons as in the previous 
subsection. 

Now let us turn to the brackets of Z and X. Here we again have, of course, 
the same expression (23) for all rows except the last one where there must be 
(leaving aside £2) 



{a?U} 



z\ - z 2 

This leads us to the following result: 
{Zi,X 2 } = — Z\X 2 \ — X 2 Z\ 2 + 



-SMf- 1 -"! V 



a 12 , a 2 



(28) 



+ 
+ 



Z\ 



l — £ Z 2 B S ® X 'rsU^j + & j) X' + 

^— Z ^ Z rs (^Z ® ID')) + 



zi - z 2 



(/<8>/jV-l)&Zi2<f 2 



(29) 
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where we again needed a new notation - this time for the matrix 



J = 



V i 



The term containing J cancels with the first term of (28) while the combi- 
nation of the remaining terms of the LHS and Z\X<i\ gives exactly the last term 
of (29). 

We now want to calculate {det Z\ , det X 2 } using the identity (20) again. We 
have two yet unknown contributions to it which we need to compute now. The 
first one is coming from X 2 i and the second one from the last term of (29). For 
the latter, we are not going to compute it at all but just rewrite it using r\ 2 
defined above, then its contribution to the brackets is 



det Z 1 det A 2 tr 2 (I ® f N -i) V2d 2 X 2 1 



(30) 



To compute the former, first note, referring to the calculation done for tviZ 2 \ 
above that 

tribal = (XX 1 + (1 ® f N ^) & 

2 Z\ — z 2 V 

with K being the matrix K defined by (24) without the last row: 



K 



( 

1 

2 

... 

\ ... 



\ 





N - 2 

0/ 



As above, the first term doesn't contribute as tr K — and the second one 
yields 

tr 2 (t ri X 21 ) AT 1 = IfLlLfi tr (/jv-i^(A') -1 ) 
To calculate the latter trace we decompose £ over the basis of the rows of X 



N-2 



and then we have 

tr (/jv-i£(A') _1 ) 



k=l 



N-2 

)N-1,N-1 

fe=l 
A 

A • • • A a ,{ - N -V A £ 



a" 



detV 
detX' 
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So combining all non vanishing contributions we obtain 

{detZi,detX 2 } = - 1 £l±^ det Z x det Y 2 + 
2 zi - z 2 

+ detZidetX 2 tr 2 (l®/Ar-i)?72C?2^ 1 
Now putting together (27) and (31) we have 

det X 2 \ lzi + z 2 

detZi, ^ = -- dctZ 1 + 

det Y 2 J 2 Z\ — z 2 

+ [tr 2 (i ® /jv-i) M^ -1 - j/^- 1 )] det ^^rj| 

which can be further simplified if we now decompose r\ in the basis of rows of 
Y': 

N-2 

fc=i 



Then: 



and 



JV-3 



r)d' = J2 ^ k+1) + HN^a'^-V + ^ 



k=i 



tr 2 (1 ® /jv_i) (v2d 2 X^ - r^Yf 1 ) = 

N-2 



fe=i 

+ MiV-l ((X'X'-^n^.n^ - (Y'Y'^n^.n^) + 

+ ^_ 2 tr 2 (J®/ JV _ 1 )o , ( JV - 1 )A:'- 1 

= ^N-2{Z'X' )n~\,N-\ 

detZ 2 

= tiN - 2 d^x~ 2 

Taking account of this we finally get 

, , , . , 1 zi + Z9 , „ / detZ(zi)\ , , , 

{*(*), A(z 2}} = -^^) + (M^^l) (31) 

Injecting this result into (25) we see that the brackets vanish for z\ = Zi, z 2 = zj 
as they are both roots of B(z) unless i = j. In this case simply taking the limit 
z — > Zj shows that we finally have 

{zi,Wj} = S lJ z i 

4.3.3 The brackets {wi, Wj} 

The final calculation we have to do is to show that 

{wi,Wj} = 
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or, in another words, that 

{A(zi),A( Zj )} = Vi,j (32) 

To calculate the brackets of A(z) we need, of course, all the brackets of X 
and Y between themselves and with each other. And for this we are again going 
to proceed as above by reusing the results obtained for {Zi,Z 2 } for all the rows 
but the last one. 

Let us start by calculating the brackets of X. We are going to have the 
following expression in the last row: 

{£1^1,62^2} = £162 [p,di + d 2 ] 

and so, in a way similar to (29), we find 

{X\,X 2 } = — X\X 2 i — X 2 X\ 2 + 

+ {I ® f N -i)&X 12 d 2 + (/iv-i ® IfaXiid! (33) 
By repeating the steps which resulted in (30), we have: 

{detXi,detX 2 } _ 1 z x + z 2 ^dcty, detYi 



det Xi det X 2 2 z x - z 2 \ det X 2 det X x , 

+ tr 2 {{I®^N-\)iid 2 X 2 x ) + 
+ tri ((/jv_i ® J^idiXf 1 ) (34) 

with 

' Ci = tr 1 e 2 (Xf 1 l 12 ) 



+ 



C 2 = tra^i^^ai) 
In exactly the same way we also obtain 

{det Xi , det Y 2 } 1 z x + z 2 det y 



det Xi det y 2 2 z\ - z 2 det X 1 

+ tr 2 ((/«./ A r_ 1 )£ 2 y 2 - 1 ) + 

+ tn ((/jv_i ® /J/sidiXf 1 ) (35) 

with 

«! = tr 2 (d^yf 1 ) 

and a symmetrical expression for {dety,detX 2 }. 
Finally, we also find 

{det Yi,dety 2 } ,, T „ . T , n 

det det = tr 2 ( ( 7^/ JV _ l)K2 y 2 - 1 ) + 

+ ^((/w-i^/Wyf 1 ) (36) 
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Combining (34), (35) and (36) together we may calculate 



{ Ax , A 2 } 1 z\ + z 2 ( det Y x dot Y 2 det Y x dct Y 2 \ 

AiA 2 ~2zi - z 2 \dctXx + dctX 2 ~ dctXi ~ detX 2 / + 

+ tr 2 (J ® jW-i) (C 2 rf 2 X 2 - 1 - C2^ X - /S2d2^ 2 _1 + «2^ 1 ) 

+ tn(/jv_i ® /) (CidiXf 1 - CiYf 1 - mdiXr 1 + mY^ 1 ) 

where the non zero terms can be rewritten as 

triC/W-i ® I) (Mi^f 1 - ^Yf 1 ) + + tr 2 (/ ® /jv-i) (njtfe-X^ 1 - v 2 Y 2 1 ) 

with f = £ — k. 

All we have to do now is to decompose v over the rows of Y 

and notice that only the terms containing A^_ 2 contribute (the others are either 
zero or, for the last one, cancel) in order to finally obtain 

f detJTi detX 2 \ detX 1 detX 2 f (1) (2) \ 

feW = det yi dety 2 ^- 2 detZ 1 + AV_ 2 detZ 2 j 

which implies (32). 

This terminates our proof that Zi and Wj have the canonical brackets and 
means that we have succeeded in constructing the separated variables in a purely 
algebraic way. 

5 Reconstruction of l(z) 

We have seen that we could relatively easily construct a point on the Jacobian 
starting from the matrix elements of l(z). In this section we are going to show 
that, conversely, we can reconstruct these matrix elements from a set of g points 
lying on the Ricmann surface and that this gives us an isomorphism between 
the space of matrices and the J — X. 

The reconstruction of l(z) is more difficult because we cannot explicitly give 
the matrix elements in terms of the Zi and Wi for N > 2. However we are going 
to show that we can write a system of equations from which these elements can 
be found. 

5.1 Polynomials Xk 

Let us start with a point on the Jacobian (wi, Zi)\<i< g satisfying ipi(zi) = 
with ip as before. 

Then we clearly have 

f b(zi)<f>(zi) = 

\ {d(zi) + Wi) (f>(zi) = 
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and this implies that TV determinants obtained by dropping k — th row from the 
N x (N — 1) matrix made of b(zi) and d(zi) + Wi are null. First, let us consider 
those of them which contain the first row: for k = 1, . . . , N — 1 we have 



/ b(zi) 

(d(Zi) + Wi)! 



X k (wi,Zi) = det 



\ 



(d(zi) 
(d(zi) 



Wi)k+1 



\ (d(zi) + Wi) N -i J 

This determinant can be developed into a sum of monomials z\w\ with 
I < N — 2. Remembering from (11) that all elements of b(z) are of degree n — 2 
in z except for the last one which is of degree n — 1, it is not difficult to see that 
most of these monomials satisfy the condition 



k < (N - 1 - l)n - 2 



(37) 



The only terms of higher degree in z come from decomposing the determinant 
along the first row and taking the last element of it. Such terms will have the 



form 



(N-l-\)n-\ i 



and appear with coefficient equal to the leading coefficient 



of Zuv which is (being also the leading coefficient of det l{z)) a constant, i.e. lies 
in the center of the Poisson algebra. 

Hence, if we move these 'exceptional' terms to the RHS of the equations 

X k (wi,Zi) = 1 < i < g 

we obtain a system of linear non homogeneous equations for g coefficients of X k 
which we can solve and obtain (almost all) coefficients of X k in the form of (2) 
because the terms in denominator, by construction, satisfy the condition (37). 
In other words, we have expressed all X k in terms of holomorphic functions on 
the affine Jacobian of (wi, z{). 

Next we turn to the remaining (yet unused) determinant 



X N (wi, Zi) = det (d(zi) + w t ) 1 



<i<g 

It contains a term id" -1 which doesn't appear in any holomorphic differential 
and also, unlike in the previous case, we have here a lot of terms of the form 
z^ N ~ 1 ~ l ' >n ~ 1 w l i with rather complicated coefficients. 

However we may remark that, by definition of d(z) as lower bottom part of 
l(z), the coefficient of ^/^(N- 1 -')™- 1 i n the decomposition of Xn is the same as 
the coefficient of w '+ 1 z ( Ar - 1 -0™- 1 j n the decomposition of det (l(z) + w). And 
this coefficient is the dominant coefficient of the polynomial t^-i-i appearing 
in the characteristic polynomial (6) and so, according to (8) and the remark 
following it, is a constant. 

Hence after moving these constant terms to the RHS we can find all the other 
coefficients of Xn in terms of the holomorphic functions on the affine Jacobian 
in the same way as above. 

To summarize, starting from a point (u>i, Zi) on the Jacobian we have found 
all X k for 1 < k < N as holomorphic functions of it. 
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5.2 Reconstruction of l(z) 

To really reconstruct l{z) we now have to show that we can express its matrix 
elements in terms of X\(z, w), . . . , Xn(z, w) and also det (l(z) + w). 
It is convenient to consider & N x (N — 1) matrix L'(z) defined by 



L'{z) 



zb(z) 
d(z) 



( 



1 
1 

V o ... 



; (0) \ 

IN > 







o / 



Note that L(z) — 0(z n ^ 1 ) by construction as we have removed all the elements 
of degree n from it. We shall also denote by l[j{z) the polynomial kj(z) without 
the term of degree n (if any). 

We also introduce another N x (N — 1) matrix 



L(w, z) 



( o 
1 




V o 



o \ 






1 / 



+ 



zb{z) 
d(z) 



and N (N — 1) x (N — 1) matrices L' k (z) and Lk(z, w) which are obtained from 
L'(z) and and L(w, z) respectively by dropping the (fc+l)-th row (1 < k < N—l) 
or the first row for k — N. 

By definition of Lk(z,w) and Xk(z,w) we have 



det Lk(w, z) 



zX k (w,z) if k < N 
X N (w,z) if fc = TV 



(38) 



As we have already found X^ in terms of the function on the Jacobian, we'd 
like to find the matrix elements of Lk now in terms of Xf,. Of course, it cannot 
be done directly from (38) so we are going to start by extracting the part of the 
determinant linear in L' . 
For this we define by 



or, more explicitly: 



V k (w,z) 



L k (w,z) = V k (z,w) +L k {z) 



1 1N Z \ 



W 

z n 



1 rows 



N — k — 1 rows 



W 
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(there is no first row for k — N). 
Then 

detL fe (z,w) = detV k (z,w) (l + tr^" 1 ^, w)L' k (z)) + 0(L' 2 )) 
and we obviously have 

detV k (z,w) = (-lfl^w^z^-^ (39) 

(there is no (-1)^^ factor for k = N )- 

Let us now calculate the trace. We can see from the definition of 14 that 
that the inverse matrix has the form 



V k -\z,w) 



or, more formally: 



w 2 z- 3n /l<V 
-wz~ 2n /l^ 



-w 2 z n w 1 



\ 



—wz 

z~' 



-In 



!(-iy+l-3 w -(i+2-j) z n(i+l-j) 
^N-l-i^N-l-iz-niN-i) 
(_\\j-i-l w j-i-l z -n{j-i) 
o 

Then it is easy to calculate: 

tv(V k - 1 (z,w)L' k (z)) = 



1 < i < k- 1,2 <j < i + 1 
k<i<N-l,j = l 
k <i <N -l,j >i + l 
otherwise 



k-l i-l 



3 w -j- 1 v n j(T>. 



z n >(L' k ) i - j+1 , i + 



i=l j=Q 
N-l N-l-i 

+ E E (-i)^- n0+1) (^Wi, 

i—k j—0 

k — 2 k—1 

= \r',r ' 1 ]T (4)^+1,* + 

3=0 i=j+l 
N-l-k N-l-j 



j=0 



i—k 



Using (39) we can write down the part of det L k (w, z) linear in L'(z) explic- 
itly as 

AetV k {z, W )tv{V k \z,w)L' k {z)) = 



- ((4)22 + (4)33 + • • • + (4) fe -i,*-i)^- 3 * (iV - fe+1)r 

+ ((4) 21 + (4) 32 + ... + (4k fe -i) + 
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+ ((4) fe+ i, fe + . . . + (4)iv-i,iv-2 + {L' k ) ltN - 1 )w k - 1 z^ N - k - 1 > - 
- ((L' k ) k+2 M + ■■■ + (L' k ) N . ltN -s + {L' k ) hN -2)w k z^ N - k -^ ± 

± ...+ 

+ (-l)^-*- 1 ^)!,^- 2 (40) 
Remember that by definition of L' k we have 

{zbj + i (mod z n ) if i = 1 and k N 
lij+i (mod z") if i < fc 
(mod z n ) if z > fc 

so all matrix elements of b(z) and d(z) appear in (40) and, remembering (38), 
we see that we have almost achieved our goal of expressing the coefficients of 
lij(z) in terms of X k . 

To make it more precise, let us introduce yet another notation: define L'^ 

by 

n-l 

L'(z) = zkL ' ik) ( 41 ) 

fe=0 

The simple but important observation we can now make is that we can 
determine all coefficients of each L'^ from (40). Indeed, this equation has 
a triangular structure and we may first use its first line for all k to find the 
elements of the second row, then use the second line to find the elements of the 
third row using those of the second one and so on. 

Now let us leave aside the linear part (40) of det L k for a moment and 
consider the entire expression for it. More precisely, let us study L k defined by 

JV-l 

det L(w,z) = wPL K z ) 

L k are polynomials of z (only) . It is immediate to see from the definition of 
L(w,z) that the maximal degree of z in a monomial w p z q which can occur in 
det L k (w, z) is (N — p — l)n — 1 as there is always a minor in which all but one 
hj(z) are of maximal degree n containing w p (except for k = N and p = N — 1). 
So we see that deg l? k = (N — p — l)n — 1. 

Moreover, as the dominant coefficient of L k (z) comes from the minor men- 
tioned above it is an expression linear in elements of L'(™ _1 ) - in fact, it contains 
either some 1^) multiplied by U+n which is just 1 or i^jv (which is not an ele- 
ment of L') multiplied by if) with i ^ j + 1. 

According to the above, we may find entirely in terms of X k (and so in 

terms of the functions on the affine Jacobian) - and so the dominant coefficient of 
L^z) can be calculated from them as well. Now let us look at the next coefficient 
in it: clearly, for the reasons of homogeneity it may only contain terms linear in 
L' n ~ 2 or quadratic in L' n ~ 1 . The latter have been already calculated and the 
former can be deduced, again, from (40). Then we continue further with the 
third coefficient of L p k (z) and so on. 

In this way we can find the entire polynomial L k (z) in terms of X k and so 
all the matrix elements of L(w, z) can be expressed as holomorphic functions on 
the affine Jacobian because of (41). 
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The only remaining unknown coefficients of l{z) are lki(z) which don't ap- 
pear in L(w, z). We can find them from the characteristic polynomial. The part 
of it linear in lki(z) has the form 

l 11 (z)w N - 1 + l 21 {z)w N - 2 z n + ... + l N1 {z)z^ N ^ n 

from where we can directly determine the coefficients of ki(z) in terms of r(w, z) 
just as above. 

In conclusion, we have shown that it is possible to uniquely determine all 
coefficients of the matrix elements of l(z) satisfying (1) and the conditions on 
its degree from the g points lying on the Riemann surface defined by the same 
equation provided these points don't fall on the theta divisor. This achieves our 
proof. 



6 Calculation of the Euler Characteristic 

The goal of this section is to compute the g-Euler characteristic of the complex 
Cq defined below. It is the same one as discussed in [1] and [2] and we apply 
the same method (using the affine model constructed in the first parts of this 
paper) for its calculation. 



6.1 Grading 

Let us consider the polynomial ring A freely generated by all coefficients ap- 
pearing in the matrix elements of l(z). We can introduce a grading on this ring 
in the following natural way: first, prescribe degree N to z and degree Nn — 1 
to w (choice of N for the degree of z is different from the conventions used in 
[1] where z has the degree 1, but is more convenient here to avoid fractional 
degrees) . 

Next we remark that the grading of det Z(^r) must be equal to N(Nn — 1) 
because it is a polynomial of degree Nn — 1 in z with dominant coefficient equal 
to 1. So the grading of £21 (2)^32(2) ■ ■ - In n-i{z)Iin(z) is N(Nn — I) and hence 
we see that we must set degl i+ n(z) = Nn and deg/ijv = N(n — 1). From 
further considerations of homogeneity it is clear that the compatible grading 
must be defined by 

deg Uj (z) = Nn + i — j — 1 

Now writing 
we also see that 

de S l tf = N ( n - dc Sz kj +a)+i-j-l 
We may now easily calculate the character of the free ring A: 

chA = 




26 



The product on the first line can be written as 



-1 JV-l 



, ;/v„-i]! 



with the standard notations 

W = [1][2]-..[JV] 
[n] = 1-9" 

As for the second part of the product, we have 

N-l N-l N 



n< ) - nnnrin 

i=l j=l «=1 y 
n iV-1 i-1 

= nn n ^ 



Q 

a=l i=l k=-{N-i) 



= n 



r* — hi 



11 L/Vn + i-ll! 

i=l 

Combining both parts together we find that 

^ 11 [iVn + i- 2]! 



6.2 The quotient ring ^4 

Next object we are interested in is the ring T freely generated by the coefficients 
of the polynomials tk(z). From their defining equation (6) it follows that, for 
the reasons of homogeneity, the grading of tk(z) is nk — 1 and its dominant 
coefficient has grading N — k. Hence 

ch ^=nflT^v^ ( 42 ) 

fc=li=l H 

where the factor with k = N and i = 1 is excluded. 
As in [1] we define the ring Ao as 

Ao = A/{F x A), F x =J2ft ( j* ) 

k,i 



And the proposition 1 of [2] gives us its character as 



chA = ^ (43) 
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6.3 Cohomology groups 

As in the hyperelliptic case we can introduce the differential operators 
corresponding to the holomorphic forms on the Riemann surface which are 

u^-V" 1 , .„ f fc = 1,...,JV-1 



We make them act on the ring A by 

With our choice of grading for z and w we have 

deg D ik = - dcga> jfe = (Nn - l)(N - k) - Ni 
This grading is, of course, compatible with the action on A: 

D lk A {n) C A (n+de&Dik) 



The character of the ring V freely generated by all Dik is easily computed 
to be 

N-l nk-1 

chD = 



N-lnk-1 

II II T^vTrr (44) 



k=l i=l 

We also introduce the differential forms Tik dual to the vector fields and 
consider the linear spaces C h for k = 0, . . . , g spanned by forms 

Xi 1 ...i k dT il A • • • A dr ik (45) 

with x i± ik G A. 
The differential 

d = y^dTjkPjk 

ik 

acts from C k to C k+1 : as usual, we first apply D ik to the coefficients x il ,,, ik 
and then take exterior product with dr^. Thus we have a complex 

-> C° C l ■ ■ ■ -i C 9 

If we fix the degrees of to be opposite to degDik the degree of the 
differential d is and so its action descends to the graded complex 

-> c ° 4 cl ^ ■ ■ ■ ^ C 9 - 

where Cq are spanned by (45) with x^...^ G ^4o- 
6.4 Euler characteristic 

The formula (44) was the last ingredient we needed to compute 
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Using (43) we get 



N-l I rafc-1 \ Nn-1 

x n Kivn - i)fc] n [m fc] 2 ) n ^ ( 4g ) 

fc=l \ i=l / i=l 

Counting the number of g-numbers in the numerator and denominator we 
see that it is the same (N(Nn — 1) + N(N— l)/2) and hence (46) admits a finite 
limit when q — ► 1: 

X (C *) = (-1^-1)^5 (^2), 

x }~J fc J| (TVi-A:) 2 (TVn- 1)! (47) 

fe=l \ i=l / 

Let us now study how does the Euler characteristic depend on n: 

N N 

*i±i = nn - 

Xn £} 1 t\(nN + k-2 + i) 

N-l k 

x Y[Y[((nk + i-l)N - k) 

fe=l i=l 
AT fe 

x J|J|((nfc + i)iV-fc) 



fe=l i=l 



Hence when n — > oo we have 



^"+1 ^ e 2 E j^ 1 fe log fc+W log JV 

and so we see that y as a function of n grows faster than N 9 and, a fortiori, 
than 2 9 which should be its asymptotic growth rate if the highest cohomology 
functions had a form as simple as in [1]. 

We can also note that even for TV = 3 case the value of the Euler charac- 
teristic is already not compatible any more with the conjectures of that paper 
and so such simple structure of the highest cohomology group is really specific 
to the hyperelliptic case only. 

7 Conclusion 

In this paper we were able to find a new explicit algebraic construction of the 
separated variables essentially independent on the concrete integrable system 
and which can be extended to the case of quantum integrable models. Using 
this affine model we have found a relatively simple expression for the Euler 
characteristic of the Jacobian of a spectral curve. 

Examining the rate of growth of the Euler characteristic, we realized that 
the simple situation of the hyperelliptic case described in [1] doesn't hold any 
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more in general. This means that we cannot obtain the values of all observables 
in our integrable system by using just the functions with the simple poles of 
the form (2) but that we probably need to consider the the functions with the 
higher order poles as well. 
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